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A  NUMERICAL  COMPARISON  BiiTWUlIN  TWO  UNCONSTRA I  NED 
VARIATIONAL  FORMULATIONS 


ABSTRACT.  In  an  effort  to  relieve  the  often  cumbersome  burden  of 
meeting  the  requirements  on  the  end  conditions  and  to  unify  the  solution 
formulation  for  boundary-  and  initial -value  problems,  unconstrained 
variational  statements  have  been  introduced  in  conjunction  with  some 
approximate  methods.  In  the  case  of  a  boundary  value  problem,  it  is 
shewn  in  this  paper  that  two  different  variational  statements  can  be 
established:  one  is  arrived  at  by  the  use  of  the  Lagrange  multipliers, 
the  other  by  energy  considerations.  The  numerical  convergence  of  the 
solutions  associated  with  finite  element  schemes  using  one  of  these  two 
different  variational  statements  is  compared  with  that  of  the  other.  In 
the  case  of  an  initial  value  problem,  both  formulations  can  again  be 
established  when  the  adjoint  field  variable  and  the  adjoint  variational 
statement  are  introduced.  The  numerical  data  presented  here  indicate 
that  while  both  methods  generate  excellent  convergent  results  for  the 
boundary  value  problem,  the  method  of  stiff  springs  yields  results  which 
show  much  better  convergence  for  the  initial  value  problem  than  those 
achieved  by  Lagrange  multipliers. 

I .  INTRODUCTION.  In  conjunction  with  variational  methods  of  mathe¬ 
matical  physics,  it  is  often  burdensome  to  select  trial  functions  which 
are  required  to  satisfy  some  or  all  of  the  end  conditions  (sec,  for 
example,  reference  [1]).  Efforts  thus  have  been  made  to  relieve  such 
requirements  on  these  trial  functions.  Courant  and  Hilbert  have  pointed 
out  that  in  conjunction  with  boundary  value  problems,  this  can  always  be 
done  by  adding  extra  boundary  terms  in  the  variational  statement  [2]. 

Such  a  concept  has  been  applied  successfully  by  Wu  in  obtaining  solutions 
to  nonconservative  stability  problems  [3].  Wu  has  further  extended  the 
application  to  the  solutions  of  initial  value  problems  [4],  Simkins 
also  developed  unconstrained  variational  statements  for  initial  and 
boundary  value  problems  [5].  The  approaches  used  by  Wu  and  Simkins  are 
different  in  that  while  Wu,  after  Courant  and  Hilbert,  employed  the 
concept  of  a  very  large  constant  (very  stiff  spring  constant),  Simkins 
used  the  method  of  Lagrange  multipliers.  For  any  given  problem,  the 
variational  statements  arrived  at  by  the  two  approaches  are  different 
in  boundary  terms.  The  purpose  of  this  paper  is  to  compare  the  numer¬ 
ical  convergence  of  them  in  terms  of  some  simple,  but  specific,  examples. 
Both  boundary  and  initial  value  problems  are  considered. 


II .  UNCONSTRAIN  lit)  VARIATIONAL  STATEMENTS  FOR  A  BOUNDARY  VAl.tll 
PROBLEM.  Let  us  first  consider  the  transverse  vibrations  of  an  l.uler- 
Bernoulli  beam  under  axial  load.  The  differential  equation  in  nondimen- 
sionalired  form  can  be  written  as  [1]: 

>•”"  ♦  Qy”  ♦  XJy  *  0  (2-1) 

where  y  =  y(x)  is  the  transverse  displacement  of  the  beam,  as  a  function 
of  the  variable  x  along  the  column's  length  (0<x<l).  The  axial  force  is 
denoted  by  Q;  A  is  the  eigenvalue  and  a  prime  ('T  denotes  a  differentiation 
with  respect  to  x.  The  problem  is  not  defined  completely,  of  course, 
without  appropriate  boundary  conditions.  Consider  the  following  given 
conditions: 

y (0)  =  y'(0)  =  0  (2-2a,2b) 

y"(i)  =  y,M  (i)  ♦  Qy’d)  =  o  (2-2c , 2d) 

Eqs.  (2-1)  and  (2-2)  define  the  familiar  buckling  problem  of  an  F.uler 
column.  It  can  be  solved  by  methods  of  approximation  in  conjunction  with 
a  variational  statement. 

<5 1 o  -  0  (2-3a) 

where 

,  1 

I0(y)  =  7  /  t(y")2  -  Q(y')2  ♦  *2y2]dx  (2-3b) 

2  0 

Through  integrations-by-parts ,  Eqs.  (2-3)  leads  directly  the  following 

i  6I0  =  0 

=  /  (y""  ♦  Qy"  +  A2y)<5ydx 

0 

♦  y"(l)6y'  (1)  -  y"(0)<5y'  (0) 

-  [y"'(i)  ♦  Qy'(i)]<Sy(i)  ♦  [y*M  (0)  ♦  Qy'(0)]6y(0)  (2-4) 

Fq .  (2-4)  indicates  that  6I0  =  0  is  equivalent  to  the  differential  equa¬ 
tion  (2-1)  and  the  last  two  of  the  b.c.  Eq.  (2-2c,2d)  provided  that  the 
variations  6y(l)  and  6y(l )  are  chosen  arbitrarily  (thus  causing  their 
coefficients  to  vanish)  and  that  6y(0)  and  <3y' (0)  vanish  identically. 

Thus,  6I0  =  0  can  be  used  as  a  basis  of  approximate  solution  if  trial 
functions  are  chosen  which  identically  satisfy  (2-2a)  and  (2-2b).  Since 
(2-2a,2b)  must  be  "imposed"  they  are  called  "imposed  boundary  conditions". 
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The  choice  of  trial  functions  is  otherwise  arbitrary  and  convergence, 
when  achieved,  will  tend  'naturally'  toward  a  solution  satisfying  (2-2c) 
and  (2-Jd)  which  are  called  the  'natural  boundary  conditions'  of  the 
problem.  The  imposed  conditions  on  the  trial  functions  are  often  burden¬ 
some  in  the  process  of  obtaining  approximate  solutions  [1).  In  this 
paper,  two  different  methods  are  compared  which  remove  these  constraints 
on  the  trial  functions. 

The  first  approach  is  an  extension  of  the  method  of  the  1 agrange 
multipliers  in  classical  mechanics.  Suppose  one  desires  to  unconstrain 
the  boundary  condition  (2-2a)  y(0)  =  0.  The  modified  variational  state¬ 
ment  shall  take  the  form  of 

61 j  =  0  ( 2  -  5  a ) 

where 

Ij  =  I0  ♦  ay(0)  C-Sb) 

and  I0  in  (2-5b)  is  given  by  (2-3b).  Eqs.  (2-5)  then  become 

6lj  =  0  =  6Iq  ♦  aSy(O)  =  y(0)6a  (2-6a) 

1 

=  /  (> .  ♦  Qy"  ♦  A 2y)6ydx 

0 

♦  y"(l)<5y’(l)  -  y" (0) 6y 1  (0)  +  y(0)6a 
-  (y'"(l)  ♦  Qy'(l)]5y(l)  ♦  [>•"’  (0)  ♦  Qy*  (0)  ♦  a]6y(0)  (2-6b) 

It  is  clear  from  Eq.  (2-6b)  that  if  one  defines 

a  =  -  [>•”’  (0)  ♦  Qy'  (0)]  ( 2 - 7 a  1 

thus 

6a  =  -  [ 6y"'  (0)  +  Q6y'(0)]  (2-7b) 

equation  (2-ob)  becomes 

1 

61 L  =  0  =  /  (>•""  +  Qy"  +  X2y)  6ydx 

0 

♦  y"(l)Sy(l)  -  (y"(0)  *  Qy(0)]6y’(0)  -  y(0)6y'"  (0) 

-  [y'"  (i)  ♦  Qy'CDlfiyfD  (2-8) 

< 


I 


Thus,  with  a  given  in  (2-7a)  and  I.  in  (2-5b)  the  variational  statement 
61  =  0  is  equivalent  to  the  given  differential  equation  and  the  boundary 
conditions  (2-2a),  (2-2c)  and  (2-2d).  Only  (2-2b)  is  imposed  on  the  set 
of  trial  functions.  This  last  constraint  condition  can  also  be  removed 
by  the  same  process  used  above.  The  completely  unconstrained  variational 
statement  through  the  means  of  the  Lagrange  multipliers  is  the  following: 

61  «=  0  (2-9a) 


with 


Since  then 


i  =  \  I  Uy")2  -  Q(y')2  ♦  x2y2]dx 
1  0 

-  y(0)y’"  (0)  ♦  y'(0)y"(0) 


1 

61-0  »  /  (y"M  ♦  Qy"  ♦  X2y)  6ydx 
0 


f -  -  9b) 


♦  y"(l)«y*(l)  -  [y'"  (1)  ♦  Qy,(l)]6y(l) 

-  y(0)6y'«  (0)  ♦  y« (0)[6y"(0)  ♦  Q6y(0)]  (2-9c) 

It  is  clear  from  Eq.  (2-9c)  that  all  the  boundary  conditions  of  Eq.  (2-2) 
are  natural  if  the  variational  statement  of  (2-9a)  and  (2-9b)  is  used. 

The  second  approach  to  remove  the  imposed  conditions  may  be  referred 
to  as  "the  method  of  infinitely  stiff  springs".  The  functional  I0  in 
Eqs.  (2-3)  can  be  identified  with  the  nondimensionalized  energy  stored  in 
the  beam.  If  the  beam  is  considered  to  be  supported  by  two  springs  at 
x  =  0,  one  reacting  to  deflection  and  the  other  to  rotation,  the  energy 
stored  in  these  springs  can  be  included  in  the  total  energy  of  the  system. 
Thus  consider 


61  =  0  (2- 10a) 

where 

1  -  i0  ♦  y  +  \  k2ty'(°)]2 

1  1 

■  7  /  [(y4')2  -  Q(y*)2  ♦  x2y2]dx 

2  0 

kj[y(0)]2  ♦  I  k2[y’(0)]2  (2-10b) 
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where  kj  and  k^  are  the  nondimensionalized  spring  constants  for  deflection 
and  rotation  respectively  at  x  =  0.  Now  since 

6 1  =  0 

1 

=  /  ♦  Qy"  ♦  X2y)  6vdx 

0 

♦  y"(l)6y'(l)  -  [y” (0)  -  k2y ' (0) ] 6y ' (0) 

-  [y'"  (1)  ♦  Qy'(l)]«y(l)  ♦  ly'"  (0)  ♦  Qy’(0)  ♦  k]y(0)]6y(0)  (2-1 1 ) 

the  natural  boundary  conditions  are 

y"'  (0)  ♦  Qy'(0)  ♦  k1y(0)  =  0  ,  y"(0)  -  k2y* (0)  -  0 

(2-12a, 12b) 

y" (1 )  =  o  ,  y'"  (1)  ♦  Qy • Cl)  =  0  (2-12c,12d) 

It  is  clear  that  Eqs.  (2-12)  reduce  to  (2-2)  if  kj  and  k~>  become  infinitely 
large.  Hence,  the  variational  statement  (2-10)  can  serve  as  a  basis  of 
an  approximate  solution  formulation  for  the  problem  defined  by  Eqs.  (2-1) 
and  (2-2)  if  kj  and  ko  are  taken  to  be  very  large  compared  with  unity  in 
actual  computations. 

III.  UNCONSTRAINED  VARIATIONAL  STATEMENTS  FOR  AN  INITIAL  VALUE 
PROBLEM.  In  the  case  of  initial  value  problems,  similar  procedures  can 
be  used  to  free  the  initial  conditions  imposed  on  the  trial  functions. 
Examples  have  been  given  in  two  previous  papers  [4,5],  Since  initial 
value  problems  are  nonself  adjoint  by  nature,  adjoint  field  variables 
must  be  introduced  to  form  variational  statements  which  provide  the  basis 
for  approximate  solutions.  In  this  section  Lagrange  multiplier  formu¬ 
lations  will  be  compared  with  those  using  the  method  of  infinitely  stiff 
springs  -  each  method  being  used  to  relax  the  requirement  that  trial 
functions  satisfy  identically  the  imposed  conditions  arising  from  an 
initial  value  problem.  Forced  motions  of  a  spring-mass  system  is  used 
for  illustration.  The  differential  equation  for  such  a  system  can  be 
written  as 

y  ♦  u2y  «  f(t)  (3-1) 

where  y  =  y(t)  is  a  function  of  the  time  t  and  a  dot  (■)  denotes  differ¬ 
entiation  with  respect  to  t.  The  constant  ui2  =  k/m  where  k  is  the  spring 
constant  and  m,  the  mass.  The  initial  conditions  are: 

y(0)  =  a  ,  y (0)  =  b  (3-2a,2b) 
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So  generality  is  lost  if,  in  establishing  the  corresponding  variational 
statements,  one  considers  only  a  homogeneous  system.  Hence  wo  consider 
the  differential  equation: 

y  ♦  u2y  =  0  ( 3- 1 ' ) 


and  initial  condition 


>’ (0)  -  0  ,  y (0)  «  0 


(.l-r'a.r'b) 


The  fact  that  the  system  of  Eqs.  (3-1')  and  (3-2')  leads  to  a  trivial 
solution  only  is  not  of  concern  here. 


Let  z  -  z(t)  be  the  adjoint  field  variable.  First,  the  variational 
statement  obtained  by  the  use  of  Lagrange  multipliers  is  verified  to  be: 


where 


«I0  =  0 


1  _  1 

I  =  -  /  yzdt  ♦  U2J  yzdt 
°  0  0 


(3- 3a) 


(3- 3b) 


♦  y(i)z(U  -  y (0)£(0) 

Eqs.  (3-3)  lead  to 

i  =  0 

=  /  (y  ♦  w2y)6zdt  y (0)3z (0)  -  y(0)5:(0) 

0 

1  .. 

♦  /  (z  ♦  w2z)$ydt  -  z ( 1) 6y (1)  ♦  z(l)6y(l)  (3-4) 

0 


Et)-  (3-4)  states  that  6I0  =  0  is  equivalent  to  the  problem  of  I'qs.  (3-1') 
and  (3-2')  and  the  adjoint  problem  defined  by 

z  ♦  w2z  *=  0  (3-3) 


and 


z (1)  =  0  ,  z(l)  =  0 


(3-hu,t>b) 


In  as  much  as  the  variations  of  the  field  variable  6y,  6z,  etc.  are  quite 
arbitrary  and  6y  is  quite  independent  of  6:,  one  can  take  6y  =  0,  3y(l)  = 
0  and  6y(l)  =  0.  Hence  the  association  of  the  problem  of  (3-1')  and 
(3-2')  with  the  variational  statement  Eqs.  (3-3)  is  established. 


Now  for  the  inhomogeneous  system  of  Eqs.  (3-1)  and  (3-2),  one  may 
similarly  verify  the  corresponding  variations  statement: 


where 


61  j  «  0  (3-7a) 


1  1 

I,(y.z)  =  -  /  yzdt  ♦  f  [w2y  -  f  (t)  ]  zdt 
0  0 


♦  y(l)z(l)  -  [y(0)  -  a]z(0)  -  bz (0)  (3-7b) 


On  the  other  hand,  when  the  "infinitely  stiff  spring”  approach  is  used 
to  treat  the  homogeneous  case,  the  variational  statement  takes  the 
following  form  [4]: 


61  =  0  (3- 8a) 

where 

1  1 

I  =  -  /  yzdt  +  ii>2/  yzdt  +  ky(0)z(l)  (3-8b) 

0  0 

Eqs.  (3-8)  result  in 

61  =  0 

=  /  (y  +  u>2y)6zdt  ♦  y(0)6z (0)  ♦  [ky(0)  -  y(l)]6z(l) 

0 

1  .. 

♦  /  (z  +  w2z)6ydt  -  z(l)6y(l)  +  [kz(l)  +  z(0)]6y(0)  (3-9) 

0 

The  differential  equations  for  the  problem  and  for  the  adjoint  problem 
are  unchanged.  The  end  condition  for  the  original  and  the  adjoint  prob¬ 
lem  are 


y (0)  =  0  ,  ky(0)  -  y(l)  =  0  (3-10a,10b) 

and 

z ( 1 )  =  0  ,  kz(l)  +  z (0)  =  0  (3-1  la, 1  lb) 

respectively,  Eqs.  (3-10)  and  (3-11)  reduce  to  (3-2')  and  (3-6)  respec¬ 
tively  as  k  becomes  infinitely  large. 
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From  Eqs .  (3-8),  extension  to  a  variational  statement  is  easily  made 
for  the  inhomogeneous  case  of  Eqs.  (3-1)  and  (3-2): 

6 1 j  =  0  (3-1 2a) 


where 


1  1 

I  =  -  /  yzdt  ♦  /  [w2y  -  f(t)]zdt 
0  0 


♦  ky (0) z (1)  -  kaz(l)  -  bz(0)  (3- 12b) 

IV.  NUMERICAL  COMPARISONS.  In  this  section,  the  two  methods  for 
the  unconstraining  of  the  coordinate  (trial)  functions  described  in  the 
previous  section  will  be  compared  numerically.  The  approximate  solu¬ 
tions  are  formulated  through  the  finite  element  discretizations. 

IV. A.  Boundary  Value  Problem.  The  example  given  in  Section  II 
shall  be  usedT  The  set  of  Eqs.  (T-l)  and  (2-2)  constitute  an  eigen¬ 
value  problem.  Using  the  method  of  Lagrange  multipliers,  the  associated 
variational  statement  is  given  in  Eqs.  (2-9)  which  can  also  be  written 

as 

1 

51  =  0  =  J  (y"6y"  -  Qy'Sy'  *  A2y6y)dx 
0 

-  y(0)5y"'  (0)  -  y'»  (0)6y(0) 


♦  y ' (0) 6y" (0)  ♦  y"(0)6y'(0)  (4-1) 

In  applying  the  standard  finite  element  discretization  the  beam  is 
divided  into  K  equal  elements.  Denoting  the  local  coordinate  by  C . 
one  has,  for  the  m-th  element: 

t  =  C(m)  =  Kx  -  m  +  1  (4 -2a) 

d£  =  Kdx  (4-2b) 

Thus,  in  terms  of  local  variables,  Eq.  (4-1)  becomes 

K  1 

61  =  0  =  l  J  [KV^'V^"  -  QKyWfiy^'  +  —  y^6ylm^]dC 

m=l  0  K 

-  K3y(1)(0)6y(1)"'  (0)  -  K3yfl)"'  (0)6y(0) 

+  K3y(1) '  (0)6y(1)"(0)  4-  K3y(1)"(0)6y  W  (0)  (4-3) 


_ 
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Bj  -  a'"  (O)a(O) 


(1  0  0  0]=  [i:  0  0  0 


6  0  0  0 


-12  0  0  0 


6  0  0  0 


(4  -  8b) 


=  a"(0)a'  (0)  =  f-6l  10  1  0  0]  ■  fl)  -6  0  o] 


0-400 


0  6  0  0 


0-200 


(4- 8c) 


Now,  I.q.  (4-7  )  can  be  assembled  into  a  global  matrix  equation 

61  -  6YT [ K  ♦  X2M]Y  =  0 


(4-  9) 


where 

VT  .  (1)  m  y  (1)  (1)  Y  m  r.) . Y  w  x  (K>)  (4-m 

1234  3  4  34’ 

The  details  of  obtaining  the  global  matrices  K  and  M  have  been  given 
e 1 sewhere  m  and  will  not  be  repeated  here. 

Since  6Y  in  (4-9)  is  unconstrained,  the  equation  reduces  to 


(K  ♦  X2M)Y  =  0 

which  will  be  solved  for  the  eigenvalues  X2. 


(4  11) 


When  the  method  of  infinitely  stiff  springs  is  used,  the  variational 
statement  is  given  by  F.qs.  (2-10),  which  can  also  be  written  as 


61  =  0  =  /  (y"6y"  -  Qy'6y'  ♦  Xzy6y)dx 
0 


♦  k,y(0)6y(0)  ♦  koy'(0)6y'(0) 
1  *> 


(4- 12a) 


K  1 


-  I  /  ( KV(l)Vm),,-QKyw,*yw,^yw*yw)d{ 


m=l  0 


♦  kiy(1)(0)6y(1)(0)  ♦  k:K:y(1)' (0)6yllV (0)  (41 2b) 
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Or, 


wheie 


K  t 

51  =  0  ■  [  6Y(m)  [K3C  -  ()KB  ♦  —  A]Y(m) 
m= 1  '  ”  *  - 

T 

♦  6Y(1)  [kjB3  ♦  k,K?B4]Y(1) 


B3  =  a(0)aT(0) 


1  O'  0  0 
OOOO 
0  0  0  0 
0  0  0  0 


(4-13) 


(4- 14:0 


B.  =  a'(0)a'T(0) 


0  0  0  0 
0  10  0 
0  0  0  0 
0  0  0  0 


(4- 14b) 


As  before,  Fq.  (4-13)  can  be  assembled  into  a  global  equation 

61  =  0  =  5YT (K  A2M)Y  (4-  IS) 

so  that  the  eigenvalue  A2  can  be  solved  from 

(K  ♦  A2M)Y  =  0  (4-16) 

Numerical  data  for  the  vibration  frequencies  of  a  cantilevered 
column  are  given  in  Tables  I  and  II  for  both  the  method  of  Lagrange 
multipliers  and  the  method  of  infinitely  stiff  springs.  As  show  in 
these  Tables,  both  methods  display  excellent  convergence. 


In  the  case  of  the  stiff  spring  method.  Tables  I  and  II  also 
indicate  that  the  greater  values  of  kj  and  k^  may  not  give  more  accurate 
results,  although  all  the  results  are  good  when  kj  and  k->  are  sufficiently 
large.  This  point  is  further  demonstrated  by  the  computations  shown 
in  Table  III.  Since  greater  values  of  k,  and  k2  mean  that  the  prescribed 
end  conditions  are  more  accurately  satisfied.  Table  III  suggests  that 
forcing  the  solution  to  greater  accuracy  at  one  point  may  cause  a  decline 
in  overall  acceptability  of  the  results  as  evidenced  by  the  declining 
accuracy  of  the  eigenvalue.  This  same  conclusion  was  first  presented 
in  [1]. 
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IV.  B.  An  Initial  Value  Problem.  For  our  numerical  comparisons  in 
the  case  of  an  initial  value  problem,  we  shall  consider  the  one  defined 
by: 

D.E.:  my  ♦  ky  =  fQcosu>ft ,  0<t<T  (4-17) 

I.C.:  y (0)  =  a  ,  y  (0)  =  b  (4-18a,18b) 

The  specific  values  of  the  constants  m,  k,  fo,  u>f,  a,  b  and  T  will  be 
given  later.  The  upper  limit  of  the  time  interval  T  can  take  any 
positive  value  other  than  infinity.  Before  one  applies  the  variational 
formulation  given  in  Section  III,  it  will  be  convenient  to  normalize  the 
time  variable  t  with  respect  to  T.  Thus  let 

t  *=  t/T  ,  t  ■=  Tt,  dt  =  Tdt 

y(t)  •  y(T)  ,  *  .  ,  iiz.  -L^i 

dt  T  dx  dt2  T2  dt2 

Also  define 

f  T2 
-  o 

u>  =  10T  ,  f  =  - 


wf  =  “f^  *  a  =  a ,  b  c  bT 

With  these  new  parameters,  Eqs.  (4-17)  and  (4-18)  become 

D.E.  -5-  ♦  i*j2y  =  f  cos(w>^.t)  ,  0<t^1 

I.C.  y (0)  *  a  ,  y(0)  *  b 

Now  we  are  ready  to  apply  the  formulations  given  in  Section  111.  We 
shall  first  consider  the  solution  formulation  by  the  method  of  Lagrange 
multipliers.  Comparing  Eqs.  (4-22)  and  (4-23)  with  (3-1)  and  (3-2), 
one  observes  that  the  variational  statement  follows  that  of  Eqs.  (3-7). 
Or, 


(4-19) 

(4-20) 

(4-21) 

(4-22) 
(4-23a , 23b) 


where 


«I  *  0  (4-24a) 

I  .  .  1 

1  =  -/  y  s  dt  ♦  |  (u2y  -  f  cos(u)ft)]z  dt 
0  0 
•  • 

♦  y(l)*(l)  -  y(0)z(0)  ♦  a  i(0)  -  b  z  (0)  (4-24b) 
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Since  6y  and  fiz  are  quite  independent  of  each  other,  one  can  set  Sy 
0  in  l-qs.  (4-24)  and  obtain 


1  a  a  1  .  .  .  1 

(61),  ■  "  /  y6:dx  ♦  J  w2y6zdt  -  /  f  cos (u<.x)dzdT 

6y*0  n  n  n  i 


♦  y(l)6z(l)  -  y (0) 6z  (0)  ♦  a6z (0)  -  b<5z(0)  =  0  (4 - 2S) 

The  same  process  of  finite  element  discretization  used  for  the  boundary 
value  problem  in  the  previous  subsection  can  be  employed  here.  The  same 
shape  functions  and  generalized  coordinates  are  also  used.  In  terms  of 
the  element  variables,  £,  defined  before,  except  now  that 


£  »  kt  -  m  ♦  1 


etc.,  liq .  (4-25)  becomes: 


(61)  -  0  -  l  5z(m)  [_Kj  a'a'Td£  ♦  —  /  aaTd(.]Y(m) 

6>’“°  m*l  0  *  *  K  o  -- 


(4  -  2t>) 


♦  6z(K)TKa(l)a'T(l)Y(K)  -  6z(1)  TKa' (0)aT(0)Y(1) 

-  I  6z(m)T  1  /  cos [  HI  (C  ♦  m  -  1 ) ] a  dC 

m=I  ”  K  o  * 

♦  5z(1)TaKa' (0)  -  6z(1)Tba(0) 


(4-27) 


I  8d”>T[-KB  .  £a]Y<">  .  4i‘«Tk» ,.»<«  -  Y(l) 

-  I  6z(m)T  \  F(m)  ♦  6z(1)T(aka'(0)  -  ba(0)]  =  0 
m=  1  K  - 

where  A,  B  have  been  defined  in  F.qs.  (4 -8a)  and 


(4-28) 


5  =  a(l)a,T(l)  «  Tol  (0  00  1]=  T0000 
0  0  0  0  0 


0  0  0  1 


0  0  0  0 


(4-20.0 


1 


B6  -  a'(0)a‘(0) 


~0~ 

[1000]= 

"o  0  0  o" 

(4 -29b) 

1 

10  0  0 

0 

0  0  0  0 

_0_ 

_0  0  0  0_ 

I 

J 

F^-1  *  /  cos[  -l  (d  +  m  -  l)a  d£ 
0  K 


1 


(4-29c) 


In  terms  of  global  generalized  coordinates  Y  and  Z  defined  by 

yT  s  [Y^l)  Y2CD  y^(l)  y^Cl)  y^(2)  y^(2) . y^(K)  y^(K)}  (4-30a) 


and 


ZT  -  [2l(1)  Z2<15  Z3<‘>  Z4(1)  Z3(2)  z„'2> . Z3<K>  Z4'K’]  (4-30b) 


Eq.  (4-28)  can  be  assembled  as  before  into  the  matrix  equation 

6ZT[KY  -  F]  =  0  (4-31) 

Or,  since  61  is  not  constrained  in  any  way, 

KY  =  F  (4-32) 

which  can  be  solved  for  Y. 

When  the  method  of  infinitely  stiff  springs  is  used,  the  variational 
statement  must  be  modified  according  to  Eqs.  (3-12).  Thus,  the  finite 
element  discretization  begins  with 

(61)  =  0 
6y=0 

1  .  .  1 

=  -  /  y6zdt  ♦  /  [w2y  -  F  cos(u£t)]z  dr 
0  0 


♦  ky(0)z(l)  -  kaz(l)  -  bz(0) 


(4-33) 


1  1 


4 
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Hence, 


I  «2(“)T[-K/‘  a'a'V  *  id  / 

m=l  O'*  K  o  -- 


2(K)1a(l)aT(0)Y<:1^ 


Or, 


♦  k6Z 


-  [  5Z(m)  1  /  cos[  (£  +  m  .  i) ]a  d£ 

m=l  -  k  o  K 

-  k6Z(K)Taa(l)  -  6Z(1)Tba(0) 

l  6Z^T(-KB  ♦  —  A]yfm)  ♦  6Z{K^kB7Yt1) 
m=l  K  -  -  -  -7. 


K 

l  <5Z 

m=  1  ’ 


lra)T  f  (m)  (K)t,  - 

jr  l  -  «Zl  ;  kaa(l)  -  6Z 


(I)Tba(0) 


w,1ere  A,  B,  F have  all  been  defined  before  and 


-  a(l)aT(0)  = 

~0~ 

fl  o  0  0]  = 

'0  0  0  0  ' 

0 

0  0  0  0 

1 

1 

10  0  0 

_0_ 

_0  0  0  0_ 

(4-34) 


(4-35) 


(4-36) 


Now  as  with  Eq  (4-28),  here  Eq.  (4-35)  can  be  assembled  in  a  clobal 
equation  in  the  form  of  Eqs.  (4-31)  and  (4-32)  and  be  solved.  8  1 


The  specific  problem  considered  is  as  follows: 


with 


my  ♦  ky  =  fo  cos(wft)  ,  0<t£T 


>■(0)  =  yD  and  y  (0)  =  yD 
The  numerical  values  of  the  parameters  are: 

m  =  1.0,  k  =  1.0,  fQ  =  l.o,  Wf  =  0.5 

y0  =  1>0»  yi  =  1*0 


IS 


The  plot  for  the  forcing  function  f0  cos(iapt)  and  the  exact  solution 
y(t)  is  shown  in  Figure  1.  The  numerical  solutions  of  the  problem 
using  both  the  method  of  Lagrange  multipliers  and  the  method  of  stiff 
springs  are  given  in  Tables  IV  through  IX. 

Tables  IV  through  VI  show  the  stiff  spring  method  generates 
excellent  convergent  results  for  various  lengths  of  intervals  of 
solut ion. 

The  results  using  the  method  of  Lagrange  multipliers  are  shown 
in  Tables  VII  through  IX.  Table  VII  shows  that  for  moderately  long 
intervals,  the  convergence  at  the  initial  point  is  non  existent 
although  it  improves  remarkably  away  from  the  initial  point.  This 
data  may  lead  one  to  doubt  whether  the  method  of  Lagrange  multipliers 
works  at  all  in  treating  i.v.  problems.  However,  when  the  length  of 
the  interval  of  solution  is  reduced,  as  shown  in  Tables  VIII  and  IX,  it 
is  clear  that  the  results  do  converge.  Hence,  both  methods  generate 
convergent  results.  The  length  of  interval  used  in  the  Lagrange 
multipliers  approach  is  so  small  compared  with  the  stiff  spring  method 
for  comparable  convergence  that  the  practical  value  of  the  former  is 
doubtful  in  treating  initial  value  problems  when  finite  element 
discretization  is  employed.  Simkins  [4]  has  shown,  however,  that  when 
global  approximating  functions  are  employed,  (consisting  of  higher 
ordered  polynomials),  very  good  results  can  be  achieved  over  an 
acceptable  interval  of  solution. 

V.  CONCLUSIONS.  From  the  numerical  data  presented  in  this  paper, 
the  following  conclusions  are  suggested: 

1.  Both  the  method  of  Lagrange  multipliers  and  the  method  of 
stiff  springs  generate  convergent  results. 

2.  In  the  case  of  boundary  value  problems,  both  methods  give  excellent 
results  and  equally  fast  convergence.  The  method  of  stiff  springs 
appears  to  be  easier  to  use  and  more  general  in  a  practical  sense. 

3.  For  initial  problems  discretized  by  finite  elements  (piecewise 
continuous  third  order  polynomials),  convergence  of  the  Lagrange  multi¬ 
plier  method,  as  compared  to  the  method  of  stiff  springs,  is  so  infer¬ 
ior  as  to  be  of  dubious  practical  value.  (This  statement  does  not  apply, 
however,  where  a  global  discretization  is  employed  using  higher  ordered 
(e.g.  8th  order  (4j)  polynomials  continuous  over  the  entire  domain  of 
integrat ion. ) 
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EXTERNAL  DISTRIBUTION  LIST  (CONT) 


NO.  OF 
COPIES 


COMMANDER 

US  ARMY  RESEARCH  OFFICE 
P.O.  BOX  1211 

RESEARCH  TRIANGLE  PARK,  NC  277C9 
COMMANDER 

US  .ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2 BOO  POWDER  MILL  ROAD 
ADELPHIA,  MD  207«3 

DIRECTOR 

US  -ARMY  INDUSTRIAL  BASE  ENG  ACT 

ATTN:  DRXPE-MT 

ROCK  ISLAND,  LL  61201 

CHIEF,  MATERIALS  BRANCH 
US  ARMY  R&S  GROUP,  EUR 
BOX  65,  FPO  N.Y.  09510 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 
ATTN:  CHIEF,  MAT  SCIENCE  DIV 
DAHLGREN,  VA  221&S 

DIRECrOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  DIR,  MECH  DIV 

CODE  26-27  (DOC  LIB) 
WASHINGTON,  D.C.  20375 

NASA  SCIENTIFIC  *  TECH  INFO  FAC 
P.O.  BOX  8757,  ATTN:  ACQ  BR 
BALTIMORE/WASHINGTON  INTL  AIRPORT 
MARYLAND  21240 


COMMANDER 
DEFENSE  DOCU  CEN 
1  ATTN:  DDC-TCA 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 

METALS  A  CERAMICS  INFO  CEN 

1  BATTELLE  COLUMBUS  LAB 

505  KING  AVE 

COLUMBUS,  OHIO  43201 

MPDC 

13919  W.  BAY  SHORE  DR. 

1  TRAVERSE  CITY,  MI  49684 

MATERIEL  SYSTEMS  ANALYSIS  A CTV 
ATTN:  DRXSY-MP 
1  ABERDEEN  PROVING  GROUND 
MARYLAND  21005 


1 

1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  BENET  WEAPONS  LABORATORY 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 
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12 
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EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 

RESEARCH  A  DEVELOPMENT 

ATTN:  DEP  FOR  SCI  A  TECH  1 

THE  PENTAGON 

WASHINGTON,  D.C.  20315 

COMMANDER 

US  ARMY  MAT  DEV  A  READ.  COMD 
ATTN:  DRCDE  1 

5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 

COMMANDER 
US  ARMY  ARRADCOM 

ATTN:  DRDAR-TSS  2 

DRDAR-LCA  (PLASTICS  TECH  1 
EVAL  CEN) 

DOVER,  NJ  07801 

COMMANDER 
US  ARMY  ARRCOM 

ATTN:  DRSAR-LEP— L  1 

ROCK  ISLAND  ARSENAL 
ROCK  ISLAND,  IL  61299 

DIRECTOR 

US  ARMY  3ALLISTIC  RESEARCH  LABORATORY 
ATTN:  DRDAR-TSB-S  (STINFO)  1 

ABERDEEN  PROVING  GROUND,  MD  21005 

I 

COMMANDER 

US  ARMY  ELECTRONICS  COMD 

ATTN:  TECH  LIB  1 

FT  MONMOUTH,  NJ  07703 

COMMANDER 

US  ARMY  MOBILITY  EQUIP  RAD  COMD 
ATTN:  TECH  LIB  1 

FT  HELVOIR,  VA  22060 


NO.  OF 
COPTER 

COMMANDER 

US  ARMY  TANK-AUTMV  RAD  COMD 


ATTN:  TECH  LIB  -  DRDTA-UL  1 

MAT  LAB  -  DRDTA-RK  1 

WARREN,  MICHIGAN  48 090 

COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  EN®  DEPT  1 

WEST  POINT,  NY  10996 

COMMANDER 
REDSTONE  ARSENAL 

ATTN:  DRSMI-RB  2 

DRSMI-RRS  1 

DRSMI-RSM  1 

ALABAMA  35809 


COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SARRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61202 

COMMANDER 

HQ,  US  ARMY  A7N  SCH 

ATTN:  OFC  OF  THE  LIBRARIAN  1 

FT  RUCKER,  ALABAMA  36362 

COMMANDER 

US  ARMY  FGN  SCIENCE  A  TECH  CEN 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  MATERIALS  A  MECHANICS 
RESEARCH  CENTER 

ATTN:  TECH  LIB  -DRXMR-PL  2 

WATERTOWN,  MASS  02172 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  B3NET  WEAPONS  LABORATORY, 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATER7LIET,  N.Y.  12139,  OF  ANY 
REQUIRED  CHANGES. 


